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We provide a quantum field theoretical derivation of the Abraham-Lorentz-Dirac (ALD) equation, 
describing the motion of an electric point charge sourcing an electromagnetic field, which back-reacts 
on the charge as a self-force, and the Mino-Sasaki-Tanaka-Quinn-Wald (MSTQW) equation describ- 
ing the motion of a point mass with self-force interacting with the Unearized metric perturbations 
caused by the mass off an otherwise vacuous curved background spacetime. We regularize the for- 
mally divergent self-force by smearing the direct part of the retarded Green's function and using a 
quasilocal expansion. We also derive the ALD-Langevin and the MSTQW-Langevin equations with 
a classical stochastic force accounting for the effect of the quantum fiuctuations in the field, which 
causes small fluctuations on the particle trajectory. These equations will be useful for studying the 
stochastic motion of charges and small masses under the influence of both quantum and classical 
noise sources, derived either self-consistently or put in by hand phenomenologically. We also show 
that history-dependent noise-induced drift motions could arise from such stochastic sources on the 
trajectory that could be a hidden feature of gravitational wave forms hitherto unknown. 
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I. INTRODUCTION 

Recent years have seen rapidly growing interest in the 
issues and problems of radiation reaction on particle mo- 
tion in a curved spacetime as a result of the need to 
produce accurate gravitational waveforms for the pro- 
posed space-based gravitational interferometer LISA. For 
generic particle-field interactions the back reaction of 
emitted radiation on the particle results from two types 
of interactions with the radiation. The first is a purely 
local effect and describes the force on the particle when 
it emits radiation. The second is non-local in nature re- 
sulting from previously emitted radiation back-scattering 
off the background spacetime curvature and interacting 
with the particle at a later time in some other position. 
The effects of the emitted radiation appear as a force on 
the particle, the self-force, that changes the particle's mo- 
tion from a simple geodesic in the background spacetime. 
DeWitt and Brehme 0] were the first to study the elec- 
tromagnetic self-force in a curved spacetime. The cor- 
responding equation depicting the dynamics of a point 
mass experiencing self-force from interactions with the 
(linear) metric perturbations it generates off a vacuum 
background spacetime was first obtained by Mino, Sasaki 
and Tanaka [2| and later confirmed by Quinn and Wald 
using axiomatic methods. For an excellent review of 
these derivations and other issues concerning the motion 
of particles in a curved spacetime see 0|. 

Meanwhile, the introduction and development of world 
line path integral techniques (see ^5] and references 
therein) allows one to derive the equations of motion 
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of (relativistic) particles interacting with quantum fields 
within a self-consistent formalism. Augmented by the 
paradigm of open quantum systems provides one with 
a self-consistent means for studying the effects of quan- 
tum field fluctuations on the particle's motion. One can 
show ah initio how these fiuctuations manifest as classical 
stochastic forces (see 0, SSl for the case of accelerat- 
ing detectors (atoms) and [HllflHTill for moving charges). 
This approach can be applied to the relation between 
Unruh radiation and ordinary radiation in accelerating 
detectors (see [3. lllllT^) . 

In a previous paper |l3l | we use the world line influence 
functional approach to obtain the equations of motion 
for a scalar charge coupled to a quantum scalar field in 
a curved spacetime. Our results for the classical motion 
agree with Quinn We also derive a scalar ALD- 

Langevin equation and show that the fiuctuations in the 
trajectory depend on their past history (non-Markovian) , 
implying that secular effects may be important when 
stochastic fluctuations are included in a more complete 
description of the particle's motion. 

In this paper, we use the same approach to study the 
motions of two kinds of objects. The flrst is the dynamics 
of an electric point charge moving in a quantum electro- 
magnetic field in a curved spacetime. This generalizes 
results obtained earlier |^ from motion in a fiat space- 
time to a curved one. The second is the dynamics of 
a small mass (point particle) interacting with quantized 
linear metric perturbations off a background spacetime 
(e.g. a massive black hole). In both cases we are in- 
terested in the radiation the particle emits and how its 
back reaction alters the trajectory. Our aim is to de- 
rive the equations of motion for the quantum expectation 
value of the particle's position. If the mechanisms that 
decohere the particle is strong enough then the expec- 
tation value behaves classically. In this limit we regain 
the well-known Abraham-Lorentz-Dirac (ALD) equation 
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in a curved spacetime, which agrees with the ear her 
derivations of DeWitt and Brehme jj, and the Mino- 
Sasaki-Tanaka-Quinn-Wald (MSTQW) equation, respec- 
tively Hg. 

We observe that the self-force acting on the particle is 
formally infinite as a result of considering point-particle 
interactions with the field. This is the usual ultraviolet 
divergence arising from the assumption of a point par- 
ticle. It must be regularized to obtain a finite and un- 
ambiguous self-force. We use a method introduced and 
developed in 5, 1^ that is inspired by effective field the- 
ory, which provides a natural way to quickly identify and 
renormalize this divergence. A high energy scale is intro- 
duced into the dynamics so that certain criteria on the 
magnitude of this scale and others must be met in order 
to have a dynamics consistent with our first principles 
derivations of (semi-classical) equations. This regular- 
ization procedure is not restricted to quantum fields so 
that one can apply this to the usual problem of relativis- 
tic particle and classical field dynamics. 

After the derivation of the semi-classical (noise- 
averaged) equations of motion we incorporate the effect 
of quantum field fluctuations on the motion of the par- 
ticles by showing how they behave like classical stochas- 
tic forces. We derive an ALD-Langevin equation for the 
electric charge motion and a MSTQW-Langevin equation 
for the point mass motion. These equations describe the 
actual trajectory including stochastic components aris- 
ing from the quantum field fiuctuations when the world 
line histories are sufficiently decohered. They are useful 
for studying the stochastic motion of charges and point 
masses under the infiuence of quantum noise sources in 
curved spacetime. The same form of the Langevin equa- 
tions derived here with quantum field-induced noise can 
be used for classical stochastic forces in a phenomeno- 
logical description of astrophysical sources in stochastic 
motion. 

For such classical stochastic forces we show that the 
correlations of the stochastic noise may give rise to a 
noise-induced drift force that is analogous to the particle 
drifts encountered in plasma physics. For noise derived 
from the quantum field fiuctuations this effect should 
contain some information about the quantum statistical 
state of the particle ■ study makes it possible to 
explore in more detail these novel phenomena that may 
have some observational effects in certain atomic optics 
experiments. This drift might also intercede in obtain- 
ing more accurate calculations of gravitational waveform 
templates, although there arc likely no significant quan- 
tum fiuctuation-induccd stochastic effects in astrophysi- 
cal sources. 

This paper is organized as follows. In Section II, af- 
ter a brief summary of the world line inffuence functional 
method we describe how to obtain the semiclassical and 
stochastic semiclassical particle dynamics for a relativis- 
tic particle interacting with a tensor field. We also de- 
scribe the philosophy and implementation of the effective 
field theory inspired approach for regularizing the self- 



force. In Section III we apply our results to the specific 
example of a point charge moving through a quantum 
vector field in a curved spacetime. In Section IV we 
study the motion of a point mass moving in a curved 
background spacetime with quantum linear metric per- 
turbations. In Section V we discuss a phenomenologi- 
cal description of stochastic particle dynamics by adding 
in noise by hand. This added noise is shown to gener- 
ate a noise-induced drift in the particle motion that is 
analogous to the usual plasma drifts. In Section VI we 
discuss the various scales appearing in our approxima- 
tions in order to obtain well-defined solutions that are 
free from the usual pathologies (e.g. runaway solutions, 
pre-acceleration, etc.). Section VII contains a summary 
of our results and Appendices A, B and C provide further 
details of certain calculations used in the text. 

Most of the treatment here parallels that in > which 
the reader is referred to for more details. We use units 
where c = G = 1, unless otherwise stated, and use the 
conventions of Misner, Thorne, and Wheeler with 
signature (— , -|-, -|-). We also use the notation that an 
unprimed (primed) index refers to that component of a 
tensor field or coordinate evaluated at the point x {x') or 
proper time r (r'), as appropriate (e.g., z" = z"(t')). 



II. RELATIVISTIC PARTICLE-FIELD 
DYNAMICS 

We choose the world line coordinates z'^{t) as the sys- 
tem of interest interacting with a quantum field <&yi(x) 
considered as its environment. The subscript A denotes 
possible tensor indices. We are only interested in the 
overall influence, but not the detailed behavior of the en- 
vironment so we will coarse-grain over the quantum field 
variables. For this purpose the world line infiuence func- 
tional is used here as in Q to determine the influence of 
the field on the particle's motion in a curved background 
spacetime (with metric 5^1/ ). 

A. World line influence functional 

Assume at some initial time t^, for a given coordinate 
system, that the quantum statistical state of the sys- 
tem (particle) and environment (quantum field) is de- 
scribed by a density matrix p(zi,^Ai] z'^j^Aijti)- For 
purely technical convenience it is customary to choose 
the initial density matrix to correspond to a factorized 
state of the system and environment 

Physically, this means that all of the field modes have 
been uncorrelated with the particle at time ti. We will 
also assume that the initial state of the environment pe is 
Gaussian in the initial field configurations. (In practice, 
there may be radiation present prior to the measurement 



3 



at the initial time.) Throughout this paper we will ex- 
plicitly retain the dependence on the initial time. 

We take the action for the complete (closed) system 
to consist of an action describing the free evolutions of 
the system of interest, that of a relativistic point par- 
ticle of mass Too, the environment, that of a quantum 
field (to be coarse-grained), and their mutual interaction, 
Ss+e[z,^a] ^ Ss[z] + SE[<^A]+S,nt[z,<^A]- The System 
action is given by 



(2.1) 



where w'^ is the particle's 4-velocity and r is the world 
line's parameter (not necessarily the proper time at this 
stage). The action for the linear quantum field is given 
by 



(2.2) 



where is the kinetic operator appropriate to the 

quantum field. The interaction action is linear in the 
field variable 



-5j^(x;z)$^(a;) (2.3) 



where the current density j^(a;;z) is some functional of 
the world line coordinates whose specific form will de- 
pend upon the field chosen. 

Since the quantum field may have gauge degrees of free- 
dom it is necessary to include a gauge-fixing term to the 
action. If the gauge choice is implemented through the 
constraint Gb{^a) ~ for some appropriately chosen 
function Gb (the « denotes weak equality in the sense of 
Dirac .17J) then the procedure of Faddeev and Popov 0| 
amounts to introducing the following gauge-fixing term 
to the field action 



(2.4) 



where a is some constant that can be chosen rather ar- 
bitrarily. In this paper we will be dealing with tree-level 
fields exclusively so there is no need to keep track of the 
ghost fields in the action. For derivations in the remain- 
der of this section we assume that the function G_b is 
approximately linear in the field so that the gauge-fixing 
action Sgf is quadratic. Any nonlinear term that might 
appear in we assume to be small and negligible within 
the context of the approximations used below. 

We coarse-grain the field variables by tracing them out 
from the time-evolved density matrix to get the reduced 
density matrix for the system 

Pr{zf,z'f,tf)^ I dzidz[ I Vz Vz' ps{zi,z[]ti) 



for times tf > ti. The quantity F[z^z'] is the influence 
functional 



F[z,z'] 



(2.6) 



which describes the influence of the coarse-grained envi- 
ronment on the particle histories z and z' . (The • denotes 
spacetime integration between times ti and t f and the dif- 
ference J J and semi-sum current densities are defined 
as the difference and average oi jA{x; z) and jA{x] z'), re- 
spectively.) The Hadamard Gfj^ and retarded Green's 
functions G^^ are given by 



G^s,ix,x') = {{^Alix),^B'l{x')}) 

GTB'ix,x') = te+ix,j:){[4>Ai{x),$B'i{x')]) 



(2.7) 



(2.5) 



and satisfy V^'^G^c' = and P-^^G^^*,, = g^^,, re- 
spectively, where g^^, symbolically represents the delta- 
function source appropriate for the quantum field equa- 
tions of motion. The above are quantum fields in the in- 
teraction picture ^ai and the quantum expectation val- 
ues are given by (• • • ) = Tr^; pe{- ■ • ) with respect to the 
Gaussian initial state pe of the environment. The step 
function 6-^-{x,'E,) appearing in G^^f equals one in the 
future of the point x' and zero otherwise, and S is a 
space-like hypersurface containing x' . 



B. Semiclassical particle dynamics 

The reduced density matrix (|2.5|l for the particle can 
now be written as 

Pr{zf,z'f,tf) 

dz,dz[ pVz j\z' ps(z,,z[-t,)e^''^^^^^^^^''''^ 

(2.8) 

where the coarse-grained effective action (CGEA) is de- 
fined as 

Scgea[z, z'\ 

= Ss\z\ - Ss[z'] + 3 a ■ Grei ' + ^ J A ' GjP' ■ J^, 

(2.9) 

and provides a description of the quantum-averaged par- 
ticle dynamics. The magnitude of the influence func- 
tional, which is a measure of decoherence jl9l |. decays 
exponentially fast for two largely separated particle his- 
tories as a result of the environmentally induced deco- 
herence of the particle through its interactions with the 
quantum field fluctuations 
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If these fluctuations provide a strong enough mecha- 
nism of decoherence, which depends on the magnitude 
of , then we are justified in expanding the CGEA 
about the classical worldline f (r) for which the small 
difference of the histories zt — z'^ — z'^' , resulting from 
quantum fluctuations of the particle itself, gives the dom- 
inant contribution to the path integral in the reduced 
density matrix. A stationary phase approximation of the 
resulting reduced density matrix gives the equations of 
motion for the classical worldline 



SScGEA 



(2.11) 



which, in terms of the actions (|2.1|) and (|2.3|l . are 

m„a^ =«;/[z]$:^*(z) . (2.12) 

The proper time gauge Ufj_u^ — —1 has been chosen. 
Here, is a derivative operator obtained from the func- 
tional derivative of the current density j^{x; z) and 
is the retarded propagator of the quantum field evalu- 
ated in the initial state pE- The gauge degrees of free- 
dom of the field will (usually) be projected out by 
Wfj^. The acceleration on the left side is Duf^/dr where 
Throughout, we use an overbar to denote 
classical quantities. 



C. Stochastic semiclassical particle dynamics 

Even under the assumption of strong decoherence, 
when the classical trajectory becomes well-defined, the 
quantum fiuctuations of the field can still infiuence the 
classical motion of the particle through the particle-field 
coupling Sint and manifest as classical stochastic forces 
on the particle. We show how this comes about using the 
influence functional formalism. 

We start by invoking the relation a la Feynman and 
Vernon jl^. 



where iV is a normalization factor that is independent 
of the worldline coordinates and ^a{x) is some auxiliary 
field. Now the reduced density matrix (|2.5|l becomes 

Pr{zf,z'f,tf)^Nj dzidz[ J Vz J Vz' ps{zi,z[]ti) 
X J VUPd^A]e'/^'^'"'^^'''''^^^ (2.13) 

where the stochastic effective action (SEA) is defined as 
Ssea[z, z'; U] = ^Scgea[z, z'] - • . (2.14) 



The function £,a{x) can be interpreted as a classical 
stochastic, or noise, field HI^J with an associated (Gaus- 
sian) probability distribution functional 



(2.15) 



The fact that this is Gaussian is a direct consequence of 
the quadratic field action 5£;[<f>] and the field's linear cou- 
pling to the particle current density in 5'i„t[z, With 
respect to ^^[Cyi] this implies that has zero-mean and 
its correlator is proportional to the Hadamard function 
encoding the information about the fluctuations in the 
quantum field .27 \ 



({e^(x),^^(x')})^ - nG^^(x,x') (2.16) 

where {...)^^Nj P^A (•..)• 

Expanding the SEA around the classical solution z'^ 
and evaluating the resulting reduced density matrix us- 
ing the stationary phase approximation results in the fol- 
lowing stochastic equations of motion for the world line 
fiuctuations z'^ = — z^ 



dr' 



■ S'^^ScGEA 



Sz"' 5zf^ 
where the stochastic force is 



(2.17) 



(2.18) 



This equation describes the dynamics of small pertur- 
bations z around the semiclassical solution z that origi- 
nate from the classical, stochastic manifestation of the 
quantum field fluctuations. We can obtain a stochastic 
version of if we add ^^TT^ to the left side of l|TT7jl . 



I.e. 



5'^ScGEA 



(2.19) 



=0 



The functional derivative gives the final form of the 
stochastic equations of motion in terms of the full world 
line z^ 



(2.20) 



This equation is only valid to linear order in the fluctu- 
ations z, which is the same order as the stochastic force 
77^ [z], because we are neglecting higher order quantum 
corrections. In practice, (|2.2l)|l is expanded to linear or- 
der in z and (|2.12|l is invoked to obtain the particle fluc- 
tuation dynamics. We point out that H2.20|l is a stochas- 
tic equation of motion because observables involving z 
must be computed via the stochastic correlation func- 
tions (. . Also, notice that both the deterministic and 
the stochastic components of the self-force can push the 
particle away from its mean trajectory with respect to a 
fixed background spacetime. 

The stochastic correlation functions of the force 77^ can 
be evaluated using the correlators above. Evaluating 
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these correlation functions along the classical trajectory 
z we find that the mean of the stochastic force is zero 
and the symmetric two-point function of the stochastic 
force (noise kernel) is 



\G 



(2.21) 



It follows that the noise ry^, which is generally multiplica- 
tive and colored, depends on the particle's initial condi- 
tions through the classical trajectory and on the field's 
initial conditions via . For most fields (with the no- 
table exception of the linearized metric perturbation h^iy) 
the projection operator enforces the noise kernel to 
be gauge invariant. For equal proper times t' = t the 
Hadamard function diverges so a suitable regularization 
procedure must be used in order to make sense of (|2.21|) 
near coincidence. 

The noise kernel (j2.21|l shows that the stochastic force 
r\^i is 0(ft^/^) as are the world line fluctuations z. This 
shows that the Langevin equation (|2.17|) is between the 
tree-level and the one-loop 0{K) equations of motion and 
therefore contains information about the lowest order |2^ 
quantum fluctuations of the coarse-grained environment. 
This is true even if the environment is weakly nonlin- 
ear. This is the reason why we do not need to include 
the ghost fields in these computations since the ghosts 
first appear at one-loop order and hence provide no con- 
tribution to the stochastic semi-classical motion of the 
particle. 



D. Retarded Green's function 

In order to study the self-force we will need to know the 
form of the retarded Green's function used to construct 
the retarded field (z) . Assuming strong environment- 
induced decoherence it follows that the field is sourced 
by a classical particle current j(z) allowing us to invoke 
Hadamard's ansatz, which takes the form in a curved 
background spacetime 

GTb\^,x') = UAB'{x,x')A'/^ix,x')S+ia{x,x')) 

+VAB'ix,x')e+{-aix,x')) (2.22) 

as the sum of a "direct" part (proportional to 6+) and 
a "tail" part (proportional to 9^). This ansatz assumes 
that X and x' are connected by a unique geodesic since 
otherwise the van Vleck determinant A{x, x') will diverge 
at a caustic. The tensor Uab' involves products of the bi- 
tensor {x,x') that parallel transports a tensor at x' to 
X. The tensor Vab' satisfies the homogeneous field equa- 
tions subject to characteristic data provided along the 
forward lightcone at x' (see 0| for details). This tensor 
field is quite difficult to calculate for generic spacetimes 
so throughout the remainder we consider this object only 
formally. 

The distributions 5+ and 9+ are defined as 9+{x,T^) 
multiplying 6(a) and 9{—a), respectively, where a = 



a{x,x') is Synge's world function along the (unique) 
geodesic linking x = z{t) and x' = z(r') and is half 
of the squared geodesic distance between the two points. 



E. Regularization motivated by effective field 
theory 

The presence of the (5(cr) term in the direct part of the 
retarded Green's function indicates a divergence in the 
retarded field = G^^^ ■ when the two points x and 
x' coincide under the integral. This is the usual diver- 
gence that results from considering interactions between 
a point-particle and arbitrarily high modes of aquantum 
field. We use the regularization introduced in p (in flat 
spacetime) and extended to scalar particle dynamics in 
curved spacetimes in ■ 

In this spirit we introduce a regulator A for the field so 
that for energies much lower than the associated energy 
scale i?A we expect the semiclassical (|2.12|) and stochas- 
tic semiclassical (I2.20|) dynamics to accurately describe 
the particle motion. This approach is taken in for de- 
riving the ALD equations in flat spacetime and is easily 
extended to motions in a curved spacetime since the 
ultraviolet divergence is local in origin. 

To regularize the divergence we make the replacement 

<5+(a(a;",x'")) ^ e+{-a)5^{a) 

= 9+ix,J:)e{-a)^A^e-^^"^" 

(2.23) 

in the direct part of the retarded Green's function (or 
retarded field <&^'^*). The function (5a (cr) is smooth but 
approximates 6{a) well only if A^cr ^ 1. This inequal- 
ity will not hold if a is strictly zero so we will assume 
that (7 is small and approaching zero while maintaining 
A^ significantly larger than l/cr. For time-like separated 
points, e.g., points on a particle trajectory, 9{—a) = 1. 

Since A serves to provide a minimum resolution the 
sharp step function 0_|_(x, S) should be replaced by the 
smooth function 



_(x,S)^0a(-s) = 



21/4 



r(i/4) 



when X and x' are on a particle world line. For such 
a pair of points the dominant contribution to 5+ comes 
from those points that are nearly coincident so that the 
proper time difference s = r' — r is small but still much 
larger than the resolution scale 1/A. Using the results 
from Appendix A it is easy to show that for large A 



eAi-s) ei-s) + - 



1 



sgn(s) 



2i/4r(l/4) s3A3 

The last term falls off faster than any inverse power of 
sA, which allows us to ignore the smooth effects of the 
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regulator on the step function so that we only need to 
keep the first term 0{—s) in all that follows. 

In the next sections we study the two specific examples 
of a point charge moving through a quantum vector field 
Af^ {x) and a point (infinitesimal) mass moving through a 
quantum field of (linearized) metric perturbations, both 
on a given background spacetime. 



III. ELECTROMAGNETIC FIELD 

In this section we will study the specific example of a 
point charge e moving through a quantum vector field 
with ^A^x) = -Afiix) in a curved background spacetime. 
The actions describing the dynamics of the field (envi- 
ronment) and its interaction with the particle (system) 
are 

S,^t[z,A^] ^ j d^x^f{x-z)A^{x) (3.1) 



where the current density is 

f{x;z)^e I dT^-^^=^g'^^{x,z)u'^ 



-g 



(3.2) 



and — Wf_iAi, — W^A^. We choose the Lorentz gauge 
so that the gauge-fixing action is 



S,f[A,] = -l J d^xV^G' 



(3.3) 



where G — V^^^ is the gauge-fixing function. From 
these we can construct the infiuence functional 



where 



Gt,^tQ^{x,Y.){{A^,k''\) 



(3.5) 



and obtain the motion of the particle at the semiclassical 
and stochastic semiclassical levels. 



A. Semiclassical dynamics 

Using (|3.4|l we compute the coarse-grained effective ac- 
tion and find that 

Scgea{z, z] 

= Ss[z]z ~ Ss\z'\ + 2^ ■ G^!,li • + I j7 ■ G^" ■ V ■ 

(3.6) 



From the discussion in Section fll Bl varying this action 
yields the semiclassical equation of motion for the particle 



moO^ = ew ^\z\A''^\7:) 
where the derivative operator is 



(3.7) 



Q/3 



P L~J p -^fi 

It follows from the functional derivative of the current 
density j" integrated over a suitable test function 

S 



Szf^ir) 



d^x^r{x;z)f{x)=ew^"[z]fiz) . (3.9) 



The retarded vector field is expressed in terms of the 
retarded Green's function A^* = GJ^'^, • . Using 
Hadamard's ansatz for the two points x = z{t) and 
x' — z{t') connected by a unique geodesic, the retarded 
Green's function is 

G;'f,(x,x') = gap,{x,x')A^/^x,x')S+{<j{x,x') 

+V^f3,{x,x')0+{^a{x,x')) (3.10) 

where ga^i is the bi-tensor of parallel transport and Vap' 
satisfies the homogeneous field equation 



□14/3' - K^V^p' = 



(3.11) 



subject to certain characteristic data determined on the 
forward lightcone at x' (see Q for details). Note that 
the direct part of the Green's function gives rise to a 
divergence when the two points are light-like separated. 
In particular, the self-force in (|3.7|l requires its evaluation 
along the world line so that the only contribution coming 
from the direct part arises when the two points are equal. 
In the next section we regularize this divergence using the 
effective field theory approach. 



B. Regularization 

Introduce a regulator A such that for energies lower 
than the scale of this regulator E\ an effective point- 
particle description of the particle dynamics becomes 
suitable. Making the replacement H2.23|l in the direct 
part of the retarded Green's function we find that 1)3. 7|) 
becomes, after passing the covariant derivative through 
the integral. 



a/3r 



[(r - T').p] Lyw-^' Ai/^^A + V^^u^' 



+ga,yu-''A'/^VpSA{<y)}+ew;^^Aifiz) . 

(3.12) 

The square brackets in the first term on the right side de- 
note the coincidence limit r' ^ r of the enclosed quantity 
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along the unique geodesic connecting and z" . Using 
the relation 

[{T-T'),p]^-up (3.13) 

which can be shown to be true from ljAl|) in Appendix 
IXI this first term becomes 



e'w^'^up {u^ + [Vo,y]u^) = — w^'Ra.pu^ (3.14) 

where = 5^" projects onto a direction orthog- 

onal to the velocity u'^. Use has also been made of the 
coincidence limit of Va-y' B 



1 



(3.15) 



The last term on the right side of (|3.12|) is the tail term 
of the self-force with 



dr' WpVaY{z^,z^')u'^' 



(3.16) 




FIG. 1: Time-dependence of the coefficients appearing in 
and ll^n^ . The function (/(I) has been divided through 
by A so that it can be displayed on the same plot with C(i) 
and g(2). 



and is a history-dependent (non-Markovian) and non- 
local quantity that allows for radiation to backscatter off 
the background spacetime curvature and interact with 
the particle at a different place and time from the origi- 
nal emission. Because the Green's functions in a general 
spacetime are not known the contribution of the tail term 
to the self-force is quite difficult to compute. We will keep 
the tail term at this formal level in this paper. 

We now consider the diverging behavior of the local 
contribution to the self-force stemming from the direct 
part of the retarded Green's function. The replacement 
(|2.23|l is a good approximation so long as A^cr 1 as 
a ^ 0. This last condition allows the direct part of the 
retarded Green's function to be expanded near coinci- 
dence t' T in a quasilocal expansion so that the diver- 
gent behavior can be identified and dealt with. Appendix 
^lists the quasilocal expansions of the relevant quantities 
appearing in the proper time integrand of H3.12|l . Using 
those results the mean particle dynamics becomes 



(3.17) 



where F^^* is some external force that may be present 
and is the self-force on the charge 



— 5(2)(r)u;, \z\^ 



+e<^[z]AL7(z-) + 0(A-^) 



\ tail 



(3.18) 



where the functions C(„) and are defined in Appendix 
IbI Fig. n shows their dependence on the elapsed proper 
time r = r — Ti . These coefficients describe the time- 
dependence of the system (particle) degrees of freedom as 
the field modes with w ^ A become correlated with the 



particle after the initial time (when they were assumed 
to be uncorrelated). 

The terms proportional to A~^ are irrelevant in the 
limit of infinitely large A while those proportional to A 
are relevant but divergent. From Appendix IbI it follows 
that the only divergent term above appears as a shift 
in the mass TOq- This shift Sm — — e^(7(i)(r) renormal- 
izes the value of the (bare) mass nio to its physical value 
mpii = Trio + Sm. We therefore define mph to be the phys- 
ical mass of the point particle at rest. Notice that when 
A ^ 00 that the physical mass becomes independent of 
time and for a finite but large value of A the shift in mass 
Sm tends to a numerical constant at late times. 



In the limit that A goes to infinity C(i), (7(2) 
the classical particle equations become 



1 and 



2e^ 

mphttfj, = — u 
+ew, 



dr 

f[z]A*:^\z)- 



-w;'[z]R^piz)u^ 



ipext 



(3.19) 



This is a generalization to curved spacetime of the well 
known Abraham-Lorentz-Dirac (ALD) equation describ- 
ing the self-force of the electromagnetic field on the mo- 
tion of a point electric charge 0. 

The effective theory paradigm is a natural one to use 
when regularizing point particle divergences since one 
cannot expect classical electrodynamics (or even QED) 
to provide a valid description of the particle at arbitrar- 
ily small distances. So long as the typical distance scales 
being probed by an experiment are much larger than this 
size 1 / A the particle behaves effectively like a point par- 
ticle. 
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C. Stochastic dynamics 

We now study the effects of tlie quantum field fluctu- 
ations (manifesting as classical stochastic forces) on the 
low-energy dynamics of the particle. We begin by com- 
puting the stochastic effective action 



Ssea[z, z'] = ^Scgea[z, z'\ - ^'^ 



(3.20) 



from (|2.14() . Varying Ssea around the classical trajec- 
tory to linear order in the world line fluctuations 
and performing a stationary phase approximation in 
the reduced density matrix gives the following Langevin 
equation 



(3.21) 



where the stochastic force 77^ is related to the stochastic 
field ^(x) through 

T^M ^ew^\z\Uz) = ~-ewf\z\VpUz) . (3.22) 

As before, the retarded field A]^*^ diverges and must be 
regularized. The prescription used in the previous section 
can be carried over so that the regularized ALD-Langevin 
equation is 

{mo - e^g(i){r)) = i^;-*(T) + f^[z] + f]^[z] (3.23) 



where the (regulated) self-force is 



Dtta 



+ j{3~c^i)ir))w;'[z]R^f,{z)u^ 



-fe<'^[z]A--(z) + O(A-0. (3.24) 

One must remember that these equations are only valid 
up to linear order in the fluctuations z about the mean 
world line z. Expanding the self-force in orders of the 
fluctuations using 



fM = urn + J dr 



(3.25) 



and computing the linearization of those terms in the 
ALD-Langevin equation (|3.23|) involving the covariant r 
derivatives (e.g. a^) gives the following equation for the 
dynamics of the fluctuations 



mf^,y[z]z'' + Jf_„y[z]z'' + Kf_,^[z]z'' 



\z\T 



+ 0(A~\l^) . (3.26) 



Here we use overdots to denote d/dr. We have also used 
the semi-classical ALD equation for the mean world line 
(|3.17() in this derivation. The tensor coefficients 



Kfj^ and r^j, are given in Appendix lO 
Notice that H3.26|l is a linear differential equation for z 
with a third derivative term and contains time-dependent 



coefficients that depend on the mean trajectory, which 
has non-Mar kovian behavior. While the mean motion is 
non-Markovian the stochastic fluctuations are Markovian 
in the sense that given a mean trajectory z the fluctu- 
ations do not depend on their own past history. The 
effective mass m^^ for the fluctuations is not generically 
diagonal implying that the inertia of the fluctuations is 
not necessarily the same in all directions. This feature 
is exhibited in the other three tensor coefficients, which 
suggests that the fluctuations in one direction are linked 
with the fluctuations in the other spacetime directions. 



D. Comparisions with scalar field 

The scalar field case examined in contains inter- 
esting features that are not exhibitied in the electromag- 
netic one. There the effective mass of the scalar charge 
is found to be 



ph 



dT'V{z^,Z^') 



(3.27) 



where is the renormalized mass of the particle at late 
times. Notice that is dependent on the history of the 
particle's classical motion via the non-local contribution 
from the tail term ^ J V so that the rate of change of 
the mass is 



dm 
~d^ 



dr' Vc 



(3.28) 



whereas the renormalized mass of the electric charge is 
constant at late times. This feature of scalar charges has 
been studied in detail on some cosmological background 
spacetimes in [2ll | where they find, in the examples they 
consider, that the scalar mass can evaporate completely 
or can lose some of its rest mass only to regain it during 
later cosmological evolution. The scalar mass changes 
because the self-force ~ ed^4){z) is not orthogonal to the 
particle's four-velocity. The component of the self-force 
in the direction of it" is precisely the quantity entering 
the rate of change of above. Choosing a different 
interaction for the scalar field, say ''^ e J dTu^^d^4>^ gener- 
ates a self-force orthogonal to the 4-velocity whence the 
mass remains constant. 

This behavior of the time-dependent mass in the scalar 
case becomes manifest in the qualitative description of 
the equations describing the dynamics of the world line 
fluctuations. Whereas these fluctuations in the electro- 
magnetic case (|3.26|) are Markovian, since they do not 
depend on their own past history, the fluctuations in the 
scalar case evolve in a non-Markovian way because of the 
appearance of 

- dr ' (5" u;^ [z] (£\ z^' ) -f ^i^M [^] Ka' ) ) 

in the dynamics, where = ct/i + w^^Vq. This term 
arises from linearizing the time-dependent mass m^(T) 
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around the classical trajectory z^. Notice that the in- 
tegration over the fluctuations at past times is weighted 
against the gradient of the tail so that this behavior dis- 
appears in Minkowski spacetime. If the classical motion 
is approximately geodesic then this term would seem to 
contribute little to the fluctuation dynamics unless the 
integration provides a significant cumulative (i.e. secu- 
lar) contribution at late times. This non-Markovian term 
does not appear when the charge current couples to the 
field's gradient ^ e J dru^d^cf). 



(The factor of 2 on the right side is a result of the sym- 
metrization g^^ g^^ ). In terms of the quantum two- 
point functions the Hadamard and retarded Green's func- 
tions are 

87rG'5<7'<5'(^'^') = {{Kfi{x),h^>s'{x')}) 
SnG'-fp^,s,{x,x') = ie+{x,^){[Kp{x),h^.5'{.x')]) 

in the interaction picture. 



IV. LINEARIZED METRIC PERTURBATIONS 



A. Semiclassical dynamics 



In this section we will study the motion of a small point 
mass mo moving through a linearized quantum metric 
perturbation field with ^a{x) = hf^^{x) in a curved vac- 
uum background spacetime. The field and interaction 
actions describing the particle-field dynamics are 



1 



d^x- 



where 

r^'^(x;z) = nio 



dX 



-gh,,T^^-{x-z) (4.1) 

5\x-z) g\[x,z)g-f,(x,z)u<-u^ 



d^x^ 



(4.2) 



is the particle's stress-energy tensor. It is assumed that 
the lowest order particle dynamics describes a geodesic 
on the background vacuum spacetime so that the accel- 
eration is of the order of the (infinitesimally small) mass 
rrio- (This assumption will be used repeatedly through- 
out the remainder.) We choose the Lorentz gauge for the 
trace-reversed metric perturbation = /i^^ — 
using the gauge-fixing action 



1 

'32^ 



d'^x^/^G^'G^, 



(4.3) 



where G^j = ^'^Ip.v is the gauge-fixing function. From 
these we calculate the influence functional and find 

F[z,z'] = e-(V«)T-,. 



xe 



{r/h)T-^-{SnG'^f/''yT+^, 



(4.4) 



(The factor of 8n multiplying the Green's functions is 
a result of the numerical factors appearing in Se and 
Sint-) The retarded linearized metric perturbation is 
^0^/3 ~ l^TrGJ^^^,^, ■T'^ ^ where the retarded Green's func- 
tion solves 



9?? ^ t'"^' 



7' s' S'^jx-x') 
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(4.5) 



Using (|2.9(l and l|4.4|l we calculate the coarse-grained 
effective action 



Scgea[z,z'] = Ss[z\ - Ss[z'] 



(4.6) 



If the quantum field fluctuations provide a strong enough 
mechanism for decoherence then varying the CGEA with 
respect to zt = z^ — z^ gives the semiclassical equations 
of motion for the particle world line, which are 

a,^w;^^[z]h::^{z) (4.7) 
where the derivative operator here is 



(4.8) 

It is computed from the variation of half the stress tensor 
with respect to the world line coordinates 



1 



2(5zm(t) 



d''x^gT'^^{x-z)f{x)^w^^^[z\J{z) 



where f{x) is some suitable test function. As before, 
the quantity = g^ + u^u"" projects vectors in the 
direction of the 4-acceleration, orthogonal to u". 

It is more convenient to use the retarded Green's 
function associated with the trace-reversed perturbation 
7^1, rather than the metric perturbation /i^^ itself since 
Hadamard's ansatz can be applied directly to the for- 
mer. In particular, if G]^^^,g, is the Green's function for 
the trace-reversed perturbation then Hadamard's ansatz 
(see 3) 

gives 

G::P'''ix,x') = 2g^/g/A'/'ix,x')S+ia{x,x')) 

+V^;'''{x,x')9+{-a{x,x')) (4.9) 
and Va^j'S' satisfies the homogeneous field equation 

DVal3YS' + 2i?/^"F^.7'5' = (4.10) 
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with characteristic data provided by its values on the 
forward hght cone from x' 0. 
In terms of 7^,^, (j4.7|l becomes 



{lTp^)-\9c.^l^''iz)\^ (4.11) 



where 7'''^* — g^'^^]^^- As with the electromagnetic 
Green's function, the direct part of GJ^^^,^, gives rise to a 
divergence when the two points are light-like separated. 
The self-force in (|4.11|l requires the retarded Green's 
function to be evaluated along the particle's trajectory, 
which is time-like, so that the only contribution to the 
self-force occurs at coincidence when the two points are 
equal. In the next section, this divergence is regulated 
using our approach inspired by effective field theory. 



B. Regularization 

Introduce a regulator A through the replacement 1 
such that for particle energies much lower than this scale 
an effective description of the particle dynamics can be 
given without recourse to information about the high en- 
ergy physics that is being ignored. We will begin by ex- 
amining the first term on the right side of H4.11|l . After 
passing the derivative through the integral and remem- 
bering the assumption that a^^ = 0(mo) it follows that 



7',-i<5' 



dr' {5A(a)V^Ai/2 + Ai/2v^<5A(a)} 



(4.12) 



The first term in the square brackets on the right side 
vanishes since [(r — t^);^] = —u^ and w^^'^UaUpUf — 
0. The coincidence limit of Vafsys' is proportional to 
Rj'{ai3)S' and vanishes when contracted against three ve- 
locity factors. The last term on the right side is the tail 
term of the self force with 



(4.13) 



ji%{z) = mo / dr' W^Vc.pys'U-' u' 



describing the history-dependent effects of radiation 
emitted in the past that gets back-scattered onto the par- 
ticle by the background curvature. It is not our intent 
here to calculate this contribution to the self-force. In- 
stead, we focus on regularizing the divergence stemming 
from the direct part of the retarded Green's function. 

Recall that the replacement 5[a) (5a (cr) of (|2.23|) is a 
good approximation if A^ct ^ 1 as ct becomes vanishingly 
small, which allows the direct part of the Green's function 
to be expanded near coincidence as r' — > r in order to 
identify the divergence. Using the quasi-local expansions 



of the van-Vleck determinant and the other quantities 
appearing on the right side of 14.12|l . listed in Appendix 
IXI and recognizing that = O(too), it follows that 

^^"^l^pi^) = -™o<[^]i?aM^)"''c(i)(r) 

+w^"P^\z]^lf^{z) + 0{K-^). (4.14) 

The background spacetime is vacuous so that the Ricci 
tensor vanishes and the only contribution comes from the 
tail term 



w,, 



,"^-*(f) = <^^[f]7*1i(^) + O(A-i). (4.15) 



Performing a similar computation for wf'^{gai3j^'^^) re- 



sults in 

(4.16) 

Putting into (|4.11|) these last two equations together 
give the MSTQW equation describing the self-force on 
the point mass 



od^ = mow^f'-'iz] hifJz) + O(A-i) (4.17) 



where 



utail ^jtail ^ ^6 entail 

"a/37 ^ 7a/37 " QapQ ISe-y ■ 



(4.18) 



Notice that the self-force of the particle is determined 
solely by the tail term of the linearized perturbation h^^. 
This makes solving for the particle's motion very diffi- 
cult as its dynamics are highly nonlinear and history- 
dependent, requiring an integration over past times and 
positions. 



C. Stochastic dynamics 

We now turn to study the effects of the quantum field 
fluctuations of the linearized metric perturbations on the 
low-energy particle dynamics. In particular, we are in- 
terested in the attributes and consequences of the quan- 
tum field fluctuations manifesting as classical stochastic 
forces on the particle. From the influence functional 14.4|l 
and the discussion in Section III Gl one can show that the 
stochastic effective action for this system and environ- 
ment is 

Ssea[z, a = ^Scgea[z. z'\ ~ C^"' • J- (4.19) 

where the current density is ~ and the CGEA 

is given by H4.6|l . Assuming sufficiently strong decoher- 
ence, expanding the SEA about the mean trajectory and 
performing a stationary phase approximation in the re- 
sulting reduced density matrix gives rise to the stochastic 
particle dynamics 



(4.20) 
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where the stochastic force r?^ is related to the stochastic 
field £,^,y{x) through 

As before, the retarded field /i^'^ diverges and must be 
regularized. The prescription used in the previous section 
can be applied here so that the regularized MSTQW- 
Langevin equation reads 

TOoflp = (4.22) 

where the (regulated) self-force is 

f,[z] = m,, hi%{z) + O(A-i) (4.23) 

and the tail term is given in terms of the trace-reversed 
perturbations in H4.18|l . We do not include an external 
force here as we did for the electromagnetic case since the 
lowest order motion of astrophysical bodies is geodesic on 
the background spacetime. Also, these equations are only 
valid up to linear order in the fluctuations z about the 
mean world line z. 

Expanding the self-force in orders of the fluctu- 
ations using (|3.25() and computing the linearization of 
those terms involving the covariant r derivatives (e.g. a^) 
gives the equation for the dynamics of the fluctuations 

mf_,^[z]z'' + jf^^[z]z'' + Ki_,^[z]z'' ^r]f^[z] -f 0(A~\5^). 

(4.24) 

We have also used the semi-classical dynamics for the 
mean world line H4.17|l in this computation. The tensor 
coefficients m^jy, 7^,^ and k^i, are given in Appendix lUl In 
particular we note that m^j/ oc g^^^,. In flat spacetime one 
can show that (|4.24|) becomes in Cartesian coordinates 

moS^ = vM+0{A-\S^) (4.25) 

since the tail term vanishes identically in Minkowski 
spacetime. At 0(A°) there is no dissipation term ap- 
pearing here, which implies that the two-point function 
of z and z could grow unbounded in time in the strict 
point-particle limit A — > 00. However, if A is large but 
finite then a simple scaling argument implies that dis- 
sipation effects from the neglected 0{A~^) terms could 
begin to appear on a time scale ~ A. Dissipation from 
higher order terms arising from the nonlinearities of the 
full metric perturbation field equations might begin to 
appear on a time scale ~ 1/mo. 

As with the electromagnetic case earlier, (|4.24|l is a lin- 
ear differential equation for z. The important difference 
is that the third derivative of z gives no contribution at 
this order. Hence, only the initial position and velocity 
of the fluctuations are sufficient to obtain a unique so- 
lution. This is unlike the electromagnetic case discussed 
earlier since one needs to introduce an external force to 
obtain unambiguous, runaway-free solutions. 

The world line fluctuations may be computed from a 
complete knowledge of the mean trajectory z. While the 



semi-classical motion of z is non-Markovian generically, 
depending on its past motion to determine its present 
state, the fluctuations are Markovian in the sense that 
there is no integration over the past behavior of z. 

V. PHENOMENOLOGICAL NOISE AND 
NOISE-INDUCED DRIFT 

The noise ri^[z] in the Langevin equations describing 
the stochastic motion of the particle in (|3.2t)|) and (|4.24|) 
are obtained from coarse-graining the environment. Our 
derivations of the equations of motion for the mean and 
fluctuation world lines assume a closed system to begin 
with. This treatment has the distinct advantage that 
it can preserve the self-consistency between the system 
and the environment in considering the effects of back- 
reaction. However, in many circumstances the stochastic 
dynamics of a system is treated phenomenologically with 
a noise put in by hand. Quite generally, for a tensor 
field ^a{x) this description follows from the stochastic 
equation 

ma^ = F;-* + ^[z]+77+ (5.1) 

where F^^* is some external force and is the self-force 
on the particle arising from its (non-local and history- 
dependent) interaction with the retarded field <i>^'^*. We 
denote the phenomenological noise 77^ put in by hand by 
a + superscript to distinguish them from those derived 
from first-principles considerations. This add-on stochas- 
tic force could have a classical origin (e.g. high temper- 
ature thermal fluctuations of a bath) or it could have no 
known single identifiable origin. Furthermore, since the 
is not derived from an initially closed system it is 
likely to be inconsistent with the dynamics of the trajec- 
tory by failing to satisfy a fluctuation-dissipation relation 
for world line dis plac ements around an equilibrium tra- 
jectory (see, e.g., ^2^). In a phenomenological treatment 
one also needs to specify the noise kernel (77+ (t)?7^, (t'))^+ 
befitting the model, rather than deriving it. 

Keeping this defect in mind (i.e., no guarantee of con- 
sistency) the analysis of Sections IIII CI and IIV CI carry 
over for a phenomenological noise. Given any kind of 
noise the equations of motion for the fluctuations around 
the mean trajectory of the particle moving through a 
(classical) field subjected to the self-force from radiation 
reaction is given by 

mf,^[z]z'' + -ff^^[z]z'' + Kf_,^[z]z'' = rf_,^[z]z'' + r]'^ 

+0(A-i,z2)(5.2) 

where the tensor coefficients for the respective electro- 
magnetic and gravitational (with r^jy = 0(A~^)) cases 
are given in Appendix [Q 

It is interesting to observe that when some source of 
noise acts as a stochastic force on the particle that a 
noise-induced drifting motion results. In particular, the 
stochastic force on the particle causes it to undergo rapid 
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motions that enables the particle to experience different 
values of the inhomogeneous external fields. Averaging 
over the stochastic fluctuations results in a noise-induced 
drift that depends on the correlations of the stochastic 
force and the gradients of the external fields. Clearly, the 
noise- induced drift is a second-order effect in terms of the 
world line fiuctuations. In this discussion we do not need 
to worry about quantum corrections from higher-order 
loops in the effective action, as the noise is not necessarily 
of a quantum origin, so we can expand 13.23|l and (|4.22|) 
beyond the linear order in the fiuctuations of the particle 
trajectory. A similar behavior is encountered frequently 
in plasma physics where the time-averaged Larmor mo- 
tion results in a net velocity drift if the charge is moving 
through an inhomogeneous external magnetic field 23]. 

In order to highlight the essential physics of the noise- 
induced drift we consider the non-relativistic motion of 
an electrically charged particle and a point mass moving 
through a fiat background spacetime. Doing so allows us 
to focus on this particular issue rather than on more com- 
plex and subtle issues that arise in the fully relativistic 
problem or in curved spacetime. Furthermore, the drifts 
commonly encountered in plasma physics are described 
by charges moving non-relativistically. 

A. Electromagnetic noise-induced drift 

To find the noise-induced drift in the electromagnetic 
case, we begin with the non-relativistic limit of (|5.1|) de- 
scribing the motion of an electric point charge moving 
through a fiat background spacetime and coupled to a 
phenomenological stochastic force 

moj = e iE^{z) + eijkUjBk{z)) + fi[z] + ri'^^iz) (5.3) 

where the external force is taken to be the usual 

Lorentz force. (The Latin indices take on the values from 
1 to 3 and etjk is the Levi-Civita totally antisymmetric 
tensor with ei23 = 1.) The self-force fi — jic" hi 
contains no tail term here since there is no background 
curvature to back-scatter the emitted radiation. Substi- 
tuting z = zq -l- in (|5.3I) and expanding in powers of 
the fiuctuations gives 

m (aoi -I- (5ai) = eSi(zo) + eSz'^daEi{zQ) 
+ ^5z°-5z^dadbEi{zQ) + etijkUojBk{za) 
+ ee^jkSujBkizo) + eeijkUojSz'^daBkizo) 
+ eeijkSujSz''daBk{zo) + ^€ijkUojSz''Sz''dadtBk{zo) 

3?'^°' 3^^^' ~^ 'nti^o) + <5z"aa?7+(zo) H 

(5.4) 

We assume that the variations of the fields occur over 
distances much larger than \5z\. The world line fiuctu- 
ations are assumed very fast compared to the averaged 



motion so we expect that \5u\ ^ |u| and similarly for the 
accelerations. This allows us to make the approximation 

mSai - eeijkSu^ B'^{zq) « ?7i^(zo) (5.5) 

where the stochastic force rj^ drives the world line fiuc- 
tuations. We assume that the typical time scale of the 
fluctuations At is much larger than the time for light to 
cross the "classical" size of the particle ~ 2e^/3mc^ so 
that the radiation reaction term, which is proportional to 
Sa, can be neglected. The equation of motion for zq then 
becomes, after taking the stochastic expectation value of 
(ICTIl . 

/ \ 2e^ 

maoi « e yEi{zo) + eijkulB''{zQ)j + -^ao^ 

+ ^{Sz'^Sz') dadbE,{zo) + ee,,k{Su^z-) daB^z^) 

+ ^t,,kui{5z'^6z'') dadbB^iza) + {6z-da4{zo)) 

(5.6) 

The terms involving the stochastic averages are defined 
as the noise-induced drift force so that 

Ff"f'^ ^{Sz-Sz'') dadbE,{zo) + ee,,k{Su^Sz-) OaB^z^) 

+ ^e,jkui{5z''5z'') dadbB\zo) + ((5z'^<9a7?+(zo)) 

(5.7) 

and 

^ 2e^ • ^ 

mao eE{zQ) + e x B{zq) + + F'^"^^ [zq] (5.8) 

The first three terms of the drift are a result of the varia- 
tion of the external fields with position and the curvature 
of the external electromagnetic field lines. The last term 
of the drift results from the worldline fiuctuation away 
from zq coupling to the variation of the stochastic force. 
If the stochastic force is independent of position then this 
term will vanish identically and any contribution to the 
noise-induced drift will result from variations in the ap- 
plied electric and magnetic fields. 

To solve (|5.8|) one needs the solution to (|5.5() for the 
fluctuations in terms of the stochastic force rf"*". The 
solution is 

5u, - - (if-^)," / dt' KabV^zo) (5.9) 

where the homogeneous solution is ignored since we are 
interested in the effect of the noise on the averaged par- 
ticle motion. The integrating factor Kab is 

Kab = exp(^-^e,b™^ dt'B"(z^')^ 

= 5ab--eabm / dt'B'^{z'^) 

+ Tf^earme\„ f dt' f dt" B"^ {zt) B^ {z',' ) 
^'Ti Jti Jti 

+ ••• (5.10) 
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Integrating 5u over time gives the world line fluctuations 
dz. The terms in 6z not involving the stochastic force are 
ignored so that one is left with 

Jti Jti Jti Jti 

(5.11) 

where the distributions Fimn and Fimnp are given by 

Fimn = '^^^''^b (^Fi{zo) + eijkU^QB''{zo)j 

+ ^ei,udaB\za)5{t - h){K^y^K2rm{K-TK^sn 
F^^np = -S{t - h)5{h - t2)5^^{K^^):Ki,^ (5.12) 

where a subscript on the integrating factor K refers to 
the designated time, e.g. Ki = K{ti), etc. This ex- 
pression for the drift force is then used to solve for the 
world line coordinates zq in (|5.8|l . This is a difficult task 
given the nonlinear and non-Markovian behavior of the 
dynamics. The history-dependent contribution coming 
from the drift force requires a knowledge of zq and the 
stochastic correlation function for all times in the past. 

If the applied fields vary over a distance much larger 
than the Larmor radius then the usual drifts that occur 
in plasma physics can still be deduced from (|5.8|) . These 
drifts have been lumped into determining the motion of 
zq in order to isolate the new noise- induced drift F'^"f* 
from the usual plasma physics drifts (e.g. grad-B drift, 
etc.). 

A similar expression for the noise-induced drift deter- 
mined above can be derived in our open quantum system 
framework. The quantum noise-induced drift would cap- 
ture some of the information contained in the quantum 
fluctuations of the world line intrinsic to the particle and 
could provide information about the quantum statistical 
state of the environment and its influence on the sys- 
tem variables. However, these considerations require a 
much more delicate analysis that goes beyond the scope 
of this paper and the stochastic semiclassical approxima- 
tion adopted within. 



B. Gravitational noise-induced drift 

A similar noise-induced drift seems to exist for the 
motion of a point mass moving through linearized met- 
ric perturbations in the presence of a source of noise. 
As before we consider the particle to be moving non- 
relativistically in Minkowski spacetime and the noise to 
be phenomenological and dependent on the particle's tra- 
jectory. Proceeding along a line similar to the electro- 
magnetic case above we find that substituting z = zq + 6z 



into H5.1|l and expanding in orders of the world line fluc- 
tuations gives two equations. The flrst describes the par- 
ticle fluctuations 

mSai « riti^o) (5.13) 

and the second describes the stochastically-averaged mo- 
tion of the particle 

mao.^Ff"f'[zo] (5.14) 
where the noise-induced drift is 

F/"^*= (fc"a„,y+(zo)),, 

d„,{vti4)v^"'i4n),+ (5.15) 

The drift force then takes the same form as the electro- 
magnetic drift with no applied fields acting on the parti- 
cle. In particular, when the noise is independent of the 
particle's position there is no noise-induced drift. 

VI. VALIDITY OF THE QUASILOCAL 
EXPANSION AND ORDER-REDUCTION 

In the first half of this paper we used the influence func- 
tional formalism and various approximations to obtain 
the low-energy effective dynamics of the particle, both 
for its semiclassical (mean) and stochastic semiclassical 
(mean and fluctuations) motion. Here, the domain of 
validity of the quasilocal expansion and this semiclassi- 
cal treatment will be discussed and compared with the 
relevant scales for weak and strong radiation damping. 

In the effective field theory paradigm a regulator A 
is introduced for controlling the ultraviolet divergences 
appearing in the direct part of the retarded fields such 
that A^cr ^ 1 with a small and approaching zero. Af- 
ter expanding a near coincidence (see Appendix IXjl the 
time scale of the quasi-local expansion At = s is gov- 
erned by At ^ A^^. Recall that for elapsed times larger 
than A~^ the time-dependent coefficients C(„) and g(^n) 
rapidly approach constant values. However, at the initial 
time these coefficients vanish giving the false impression 
that the usual problems of pre-acceleration and run-away 
solutions have been resolved. This time-dependence can 
be viewed as a consequence of choosing a factorized ini- 
tial state in our derivations since it takes a time of order 
A~^ for the field modes to fully interact with the par- 
ticle as they were uncorrelated at the initial time. To 
better study the pre-acceleration and run-away solutions 
we should choose a more physical initial state than the 
factorizable one used in most treatments. 

The mean and stochastic particle dynamics are ob- 
tained here by using the Gaussian approximation to 
compute the reduced density matrix, which amounts to 
working at the tree-level in both the particle and the 
field sectors. This implies that At should be much 
longer than the time scale for creating particle pairs. 
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Ar ^ h/m = Ac where Ac is the particle's Compton 
wavelength. 

Another relevant scale appears when trying to find 
unique, physical solutions to the ALD equation H3.17|l . 
which contains a term with a third derivative of the par- 
ticle's position. As is well known, this term is responsible 
for the problematic existence of pre-accelerated, acausal 
and runaway solutions. These kinds of solutions can be 
eliminated if the self-force is weak compared to other 
forces acting on the particle. In particular, an asymptotic 
expansion in powers of tq ~ 2e^/3mc^. called the Lan- 
dau approximation or order-reduction |24| . is employed 
to obtain physical solutions that require only an initial 
position and velocity. The Landau approximation con- 
verts the ALD equation (of third order) to the so-called 
Landau-Lifshitz equation (of second order). See for 
an interesting discussion of these equations and order- 
reduction. The quantity tq is often called the "classical" 
size of the charge |^ . 

Using order-reduction, the lowest order solution is 
found by simply ignoring the self-force so that the ra- 
diation damping is assumed weak. The time-scale of the 
dynamics is then determined mostly by the external force 
so that if F^^* varies on a scale X^xt then Ar ^ Xext- In 
curved spacetime the self-force will be weak if tq <^ Ar 
and the scale associated with the spacetime curvature A^, 
is much greater than tq. 

For the electron, Ac ~ 137ro » ro and one might 
choose to set A~^ ~ Ac to justify ignoring the effects 
of electron-positron pair production from appearing in 
the semi-classical particle dynamics. For an ion or a 
larger charged body, its physical size Rq will undoubt- 
edly dwarf its "classical" size ro and Compton wavelength 
Ac so that one might choose A~^ ^ Rq in order to ig- 
nore any effect resulting from the object's finite spatial 
extent and describe the object effectively as a point par- 
ticle. Our approach would need to be augmented if we 
wished to include the effects of extended charged bodies. 
One way to do this while still maintaining the effective 
field theory paradigm is to include all possible terms into 
the point-particle action Ss[z] that are consistent with 
reparametrization and diffeomorphism invariance. This 
provides a model-independent way to parametrize the 
contributions to the dynamics from the body's size. The 
couplings of these extra terms can then be determined 
by matching this effective theory to the theory describ- 
ing the body on microscopic scales. See 26] which takes 
a similar approach to construct a framework to derive 
the post-Newtonian equations describing the motions of 
neutral (spinning) extended bodies interacting gravita- 
tionally. 

For the gravitational third derivative of the par- 

ticle's position appears at order 0(A~^,mo) so one does 
not need to implement the order-reduction procedure to 
find well-behaved solutions to (|4.17() . One may be inter- 
ested in studying the motion of a small mass black hole 
moving in the spacetime of a black hole of much greater 
mass, e.g. at least for obtaining gravitational waveforms 



for the LISA experiment. In this case, their gravitational 
radii are sufficiently different to justify describing the 
smaller body as a point particle. If the smaller body 
is a neutron star then a description that could take into 
account the star's finite extent, rapid spinning and strong 
magnetic fields might be more desirable. In such a sce- 
nario one might wish to take A~^ ~ Rns where Rns is the 
size of the neutron star. In |26j| the authors used such an 
effective point particle theory to describe the motion of 
spinning, slowly moving and weakly gravitating extended 
bodies. Among other interesting results they found that 
in certain cases the extended bodies can be treated effec- 
tively and accurately as point particles up to 0{v^'^) in 
their relative velocity. 



VII. SUMMARY 

In this paper we have given a first-principles derivation 
of the electromagnetic ALD equation for the world line 
of an electrically charged point particle interacting with 
a quantum vector field as it moves in a given back- 
ground curved spacetime. Our Eq. H3.19|l for the low- 
energy effective particle dynamics agrees with the results 
obtained earlier 0. In a similar vein, we also derived 
from first principles the MSTQW equation 0, Q describ- 
ing the motion of a point mass moving through a vacuous 
background spacetime, including the history-dependent 
self-force experienced by the body. If the quantum fluc- 
tuations in the field strongly decohere quantum histories 
of the particle then the particle's quantum fiuctuations 
can be ignored and the semiclassical motion H3.17|l and 
(|4.17l) obtained. We use an effective field theory method 
to regulate the singular behavior of the retarded field. 
For a sufficiently fine temporal resolution but still much 
coarser than the inverse of our ultraviolet regulator A, a 
quasilocal expansion can be used to obtain the contribu- 
tions relevant to the self-force and those that are irrele- 
vant in the infinite A limit. (For large and finite values 
of A the 0(A^^) terms are dependent on the particular 
rcgularization scheme implemented as well as the details 
of the microscopic theory that it must be matched to.) 
This approach provides an intuitive and effective way to 
identify and regularize the divergences appearing in the 
self-force. 

Fluctuations in the quantum field are expected to af- 
fect the particle's motion causing it to fluctuate by an 
amount z around the mean trajectory z given by the so- 
lutions to the semi-classical equation. We derived the 
correlations of the stochastic force ?7^[z] associated with 
the fluctuations of the quantum field, and a Langevin 
equation for both the electromagnetic and gravitational 
cases. The dynamics of the fiuctuation in the trajectory 
z is Markovian but depends on the non-Markovian be- 
havior of the mean trajectory z. This is to be contrasted 
with the scalar field case ([13) where it is found that 
the fiuctuations z do depend on their own history in a 
curved spacetime for the particular particle-field coupling 
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chosen. 

Instead of the noise f?/i[z] derived here, from first- 
principles as related to the fluctuations of quantum fields, 
one can consider some other classical noise 77+ suitably 
chosen to model some stochastic source in a phenomeno- 
logical description. We can still use (|3.26|) and (|4.24|) 
with rjf^[z] replaced by 77+ to study the effect of such 
noises on the fluctuations of the particle trajectory. Since 
the origin for the noise is no longer due to a quantum 
field we need not worry about keeping up to linear or- 
der in z in these Langevin equations. Instead, in flat 
spacetime and for slowly moving particles, expanding the 
solution in terms of a slowly varying world line about 
which fast fluctuations occur z = zq + Sz and taking the 
stochastic average (expectation value) of the slow mo- 
tion dynamics we show that there is, in general, a non- 
vanishing noise-induced drift force H5.11|l . H5.15|l result- 
ing from the correlations of the stochastic force. Instead 
of using phenomenological noise, we can instead derive 
this noise-induced drift effect of the quantum field fluc- 
tuations using the influence functional. Such an effect, 
while undoubtedly small, contains information about the 
quantum statistics of the field that might have obser- 
vational consequences for experiments involving the mo- 
tion of particles in the presence of boundary-constrained 
quantum fields (e.g. Casimir and Casimir-Polder effects). 
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APPENDIX A: THE QUASILOCAL EXPANSION 

In this appendix we provide the quasilocal expansions 
of the relevant quantities used in regulating the direct 
part of the retarded fields A^^ and ft.^^ in Sections UTTI 
and II VI The relations we start with have been derived in 
Appendix A of .ISj to which the reader is referred to for 
more details. The quasilocal expansion of Synge's world 
function is 



cr"(z",z" ) = -i 



2! 



^ Til 

n=l 



dr 



s-^ Da" 



(Al) 



From the identity 2(7 = ffa^'" one can show that 

+ 0{s^) (A2) 



s 
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where = Uaa" |23 |. 

We can apply these results to compute the quasilocal 
expansion of the quantities appearing in the direct part 



of the retarded fields ^J^*^* and /i™*. Following Q the 
expansion of these quantities around a, a" =0 at x" = 
z"(r) is 

A^^^iz^^z"') 

= 1 + ^ i?a/3 a" (T^ - ^ Rc.„.^ a"a^a'< + ... , 

(A3) 

Invoking the quasilocal expansion s = r' — r and using 
(lAip and l)A2|l results in the following 

Ai/2(z", z"') = 1 + ^ Rapu'^u" + 0(s3) , 
6 



d(T\ ^ OSa 



s uDa^\ dSA 3 

+ + _J_ + 0(,) 

(A4) 

where — + w^u"^ projects vectors onto a direction 
orthogonal to the 4- velocity . 

Lastly, we will need to determine the quasilocal ex- 
pansion of gaf3'U^ and gaj'-iBu'^ when calculating the 
regularized self-force. These kinds of terms appear in the 
direct part of the retarded fields and result from inte- 
grating the direct part of the retarded Green's function 
with the current density. From Q one can deduce the 
expansions around cr, cr" = at a;" = z"{t) to be 



5a/3'W^ = -CTa/S'W^ - -^cTfj'^Rajisscr'^cr^u'^ H 



1 

2' 



(A5) 



Recalling that (Jap'U^ — Da^/dr' and using the quasilo- 
cal expansions of a a and a from ljAl|l and ljA2|l . respec- 
tively, one can show that 



gaP'u'^ = Ma + SOq + ^ — + ©(s"*) , 

z dr 

ga-i';l3U' ^ -^Ra^psu'u^ + Oi-'P-) . (A6) 
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APPENDIX B: TIME DEPENDENT 
COEFFICIENT IN REGULATED DYNAMICS 

The r-dependent coefficients (r — t — Tj) appearing 
in the equations for the regulated semiclassical dynamics 
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=(„) 



(r) 



2(«-i)/4 
' ^1/2 „,! 



1 + n r^A* 



(Bl) 



of the following relations. Expanding the world line co- 



ordinates as 



implies that the velocity is 



and the acceleration 



dr 

+o{ 



2r 



Q/3 



a/3,7 



(CI) 



as 

2(n+6)/4 



(B2) 



where 7(0, b) — T{a) — T{a, h) and T(a, b) is the in- 
complete gamma function. As A or r goes to infinity 
^{a,b) r(a)- Since C(„) and g(„-) are even functions 
of the elapsed time then we obtain the same values also 
when the initial time is in the remote past. In either 
limits both C(i) (r) and 17(2) {t) approach 1 on a time scale 
r - 1/A. 



An overbar indicates that the function is evaluated at the 
mean world line coordinate z^. One can also derive an 
expression for Da^ /dr but we do not give it here. The 
indices are to be raised and lowered with the metric 



and the inverse metric is 



(C2) 



(C3) 



APPENDIX C: TENSOR COEFFICIENTS IN 
LINEARIZED STOCHASTIC DYNAMICS 

To determine the z-dependent tensor coefficients ap- 
pearing in H3.25|) and (|4.24|) we will need to make use 



Using these relations we can determine the tensor coef- 
ficients m^i, . . . r^i, appearing in the equations of motion 
for the world line fluctuations z'^. 

The z-dependent tensor coefficients appearing in (|3.26|) 
for the electromagnetic case are 



2e 



K^iu[z\ = TOph(r) [g^aT'p^^^u'^u' a"5pa,<.] ^ "T" (3 - C(i)(r)) w^°'Rap,uu'^ - ew^^"'^ A'j^^[^ 



2e2 

^9(2){r)w^a 



-ad- 

9 '^9Pi,y 



DaTi 
~d^ 



- f3-t,u 



(C4) 
(C5) 

(C6) 
(C7) 



where the convention is used that T^a[3\ — {Tap — ?)3a)/2. In flat spacetime and in Cartesian coordinates these become 
much simplified to 



4e^ 



(C8) 
(C9) 



where the symmetrization convention used is T(a/3) — {Tafj + T|3a)/2. 
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The z-dependent tensor coefRcients appearing in H4.24|l for the gravitational case are 



■a/37,1/ 



(CIO) 
(Cll) 
(C12) 



In flat spacetime, 7^1/ = k^j^ — and (|4.24l) is simply vigZ^ = ri^(z). 

We also include corrected expressions for the tensor coefficients that appear in (5.7) of Using their notation 
for the tensor coefficients we have 



mFj^{r) — I dr'V 



a -a 



gp.v - e g(2){r)w^,a^'^^u' 



mf{r)-e^ I dr'V 



gi.oXf}uU^ - — g(2){r)u(j,g''l-j^ - — 



C(i)(^) w"Ra,^ + 2u(^g"lRapu'^ 



■g{2){r)Wf,a 



2— -a) /tail 



K^4z] 



mf{r)-e' / dr'V 



-aB- 

9 '^gpi.u 



DoTi 



rt.Az] = y 5(2) (Ow^M^ 



(C13) 

(C14) 

(C15) 
(C16) 



where the physical mass of the scalar charge is 'rrJ'^{r) — rrio — e^(7(i)(r)/2 + e^C(o(7') and 



M^,[z]r + T^,[z]r + K^,,[z]r - I dr' 



r' (^S''w^[z]V.o. + S"' w^[z]V.o 
The operator = + w^Va is independent of t'. 



r^^zi'S^ + OiA-') (C17) 
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